Peierls Dimerization with Non— Adiabatic Spin— Phonon Coupling 
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We study the magnetic properties of a frustrated Heisenberg spin chain with a dynamic spin-phonon 
interaction. By Lanczos diagonalization, preserving the full lattice dynamics, we explore the non- 
adiabatic regime with phonon frequencies comparable to the exchange coupling energy which is 
e.g. the relevant limit for the spin-Peierls compound CuGeOa . When compared to the static limit 
of an alternating spin chain the magnetic properties are strongly renormalized due to the coupled 
dynamics of spin and lattice degrees of freedom. The magnitude of the spin triplet excitation 
gap changes from a strong to a weak dimerization dependence with increasing phonon frequencies 
implying the necessity to include dynamic effects in an attempt for a quantitative description of the 
spin-Peierls state in CuGeOa. 



PACS numbers: 75.10.Jm, 75.40.Mg, 75.50.Ee 



Low dimensional quantum spin systems have attracted 
considerable attention of theorists over the decades. 
Most of the remarkable features observed in these sys- 
tems are pure quantum effects uniquely due to their low 
dimensionality. Among the current experimental and 
theoretical efforts for understanding the magnetic prop- 
erties of linear chain and ladder materials the recent dis- 
covery of a spin-Peierls (SP) transition at T$p = 14.3-K" 
in the inorganic compound CuGe03 (CGO) Q has re- 
ceived particular attention. In analogy to the Peierls in- 
stability towards dimerization in quasi one-dimensional 
metals [0 the energy of spin chains is lowered by dimeriz- 
ing into an alternating pattern of weak and strong bonds. 
The SP transition is therefore driven by the magnetic en- 
ergy gain which overcompensates the lattice deformation 
energy ||. In the SP phase of CGO the copper mo- 
ments form singlet dimers along the chains and there is 
an energy gap to spin triplet excitations. Experimen- 
tally, the SP nature of the transition and the spin gap 
have been firmly established by inelastic neutron scatter- 
ing (INS), susceptibility, X-ray, and electron-diffraction 
experiments 0. 

The Peierls transition naturally involves also the lat- 
tice degrees of freedom due to the alternating distortion 
of atomic positions. A rather common situation has been 
the clear separation of the electronic or magnetic energy 
scales from the frequencies of those phonons which cou- 
ple most effectively to the SP lattice distortion. This 
adiabatic limit is e.g. realized in trans polyaccthylcne 
p| or organic SP materials §. In contrast, in CGO 
two weakly dispersive optical phonon modes involved in 
the SP transition have been identified with frequencies 
Tiuj ~ J and Tilo w 2J, where J is the superexchange 
coupling between neighboring Cu spins §. These hard 
phonons require a very strong spin-phonon coupling for a 
SP transition to occur || . A strong magnetoelastic cou- 
pling has indeed been observed in thermal expansion and 
magnetostriction experiments 0] and inferred from large 
uniaxial pressure (pi) derivatives dJ/dpt JT3]. The SP 



physics in CGO is therefore in a non-adiabatic regime. 

Despite this unusual situation previous theoretical 
studies have commonly adopted an alternating and frus- 
trated antiferromagnetic (AF) Heisenberg spin chain 
model (ll) 



H = J 
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with a static dimerization parameter <5, thus representing 
the extreme adiabatic limit of a SP chain, i denotes the 
sites of a chain with length N and Sj are 5 = 1/2 spin 
operators, a determines the strength of a frustrating AF 
next-nearest neighbor coupling. The model (|l|) contains 
two independent mechanisms for spin gap formation. At 
5 = and for a < a c the groundstate is a spin liq- 
uid and the elementary excitations are massless spinons 
|f2| . a c = 0.241 was accurately determined by numerical 
studies ]l3| , pd| . For a > a c the groundstate is sponta- 
neously dimerized, the spectrum acquires a gap and the 
elementary excitations are massive spinons Q,^5|. On 
the other hand for any finite 6 the singlet groundstate of 
(Q) is also dimerized, but the elementary excitation is a 
massive magnon Jl^,|l6| . 

From fits to magnetic susceptibility data in the uni- 
form high temperature phase Jl7| , p~0| | and from the re- 
quirement to reproduce the experimental triplet excita- 
tion gap 

0JH tne parameter set J = 1Q0K, a = 0.36, 
and S = 0.014 was estimated for CGO within the static 
model (|l|). Yet, the amplitude of the dimerization is sub- 
stantially underestimated when compared to estimates 
from structural data in the SP phase [l9[|7|] . In this letter 
we will show that due to the non-adiabatic phonon dy- 
namics a much larger dimerization is needed to achieve 
the same magnitude of the triplet gap and thus offers 
a possible route towards a quantitatively consistent de- 
scription. 

The dimerization parameter 8 in Eq. (^) may be viewed 
to result from a mean-field treatment of an interchain 
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FIG. 1. Energies of the singlet groundstate (full symbols) 
and the lowest triplet excitation (open symbols) vs. the num- 
ber of phonons used in the numerical diagonalization. The 
different symbols indicate the results obtained by including 
phonons (a) only with q = n (o), (b) with q = and q = n 
(O), and (c) at all wavevectors (□). As in all subsequent 
figures energies are measured in units of J. 

coupling which is generated by a purely elastic mecha- 
nism ]2fj| . In fact, this coupling to neighboring dimerized 
chains provides an external potential lifting the degen- 
eracy between the two possible dimerized groundstates 
plf . Thereby, soliton and antisoliton excitations form 
spin triplet or spin singlet bound states p2[ |, the former 
is the elementary massive magnon of the spin Hamilto- 
nian @) for 6 > 0. 

In this letter we choose a different starting point and 
use a Heisenberg chain model in which the spins are ex- 
plicitly coupled to the lattice degrees of freedom main- 
taining their full quantum dynamics. Specifically we 
study the magnetic properties of the Hamiltonian pi] 

H SP = JJ2 + 9( b t + M) S 4 • S i+ i +aS r S l+2 ] 

i 

+ ?kjJ2 b t b i ( 2 ) 

i 

with the local phonon creation and annihilation operators 
bf and &i, respectively. For the optical phonon with fre- 
quency u) we thus choose a dispersionless Einstein mode 
- an appropriate choice for modeling those phonons in 
CGO which are most relevant for the Peierls distortion 
. In the non-adiabatic (huj ~ J) and intermediate cou- 
pling [g ~ J) regime the spin and the phonon dynamics 
are intimately coupled and integrating out the phonons 
in order to obtain an effective renormalizcd spin-only 
Hamiltonian is in general not possible |2^| . 

This situation thus requires a numerical approach. We 
have performed exact diagonalizations on finite chains 
up to N — 16 sites with periodic boundary conditions. 



FIG. 2. Finite size scaling behavior of the triplet exci- 
tation gap A ST for a = 0.36. The lines result from the fit 
to the scaling function (^). The (almost adiabatic) phonon 
parameters g = 0.11 and ui — 0.3 were used in order to make 
contact with previous numerical studies Q . 

Since the Hilbert space associated with the phonons is 
infinite even for a finite system, we apply a controlled 
truncation procedure retaining only basis states with at 
most M phonons p4| . M is fixed by requiring a rela- 
tive error in the groundstate energy less than 10 -7 . In 
our calculations all possible phonon modes are taken into 
account. 

Fig. |l| demonstrates for a 12 sites chain how many 
phonons have to be included to achieve convergence for 
the groundstate energies of the singlet S = and the 
triplet 5 = 1 sectors. Although the specific value of a 
is not crucial for our arguments, we keep the frustra- 
tion parameter a — 0.36 in all calculations fixed to its 
most likely value in CGO. Fig. |l| contrasts the results 
obtained if only phonons with wavevector q — n or q = n 
and q = are used to the case when all phonon modes 
are included. While it has been argued in Ref. [^5| that 
keeping the q = 7r mode only should already capture 
the dominant dimerization effect of the spin-phonon in- 
teraction, we find in Fig. [l] that e.g. the singlet-triplet 
excitation gap A ST — Eq — Eq is strongly renormalizcd 
when phonons of all q are taken into account implying 
that a spin triplet excitation is accompanied by a local 
distortion in the lattice. 

The renormalization of the singlet-triplet gap A ST is 
shown in the scaling plot in Fig. ^. It is contrasted to 
the result with a restriction to the q — n phonons only. 
Despite of the limited system sizes of up to 16 sites the 
data are well fitted by the scaling function [ p7|]Tg|l 

A ST (N) = A ST + ^e- N ? N ° . (3) 

The data clearly demonstrate that the restriction to the 
7r-modes substantially overestimates the size of the gap 
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FIG. 3. Dispersion of the elementary triplet excitation 
(left scale) for a chain with N = 16 sites for different phonon 
frequencies and spin-phonon coupling constants. Experimen- 
tal data (*, right scale) are also included from Ref. [2(1. 

for N — > oo. 

The triplet dispersion shown in Fig. ^ supports the ex- 
pectation - as inferred from the neutron scattering stud- 
ies of the phonons J?]] - that the frequencies of the rele- 
vant optical phonon modes in CGO are comparable to or 
larger than the exchange coupling J. Because for Twj < J 
the phononic character of the excitation dominates lead- 
ing to a flattening of the dispersion p8| , p9| . This, how- 
ever, is incompatible with the magnetic INS data |2(f | 
which do not show any signature for a flat dispersion of 
the triplet excitation branch. We have included in Fig. || 
experimental data on CGO from Ref. j2(| which, how- 
ever, is not intended as a fit to the data - our results 
in Fig. H are for a fixed chain length with N = 16 sites 
- but rather shows the shape of the triplet dispersion. 
We stress that the upturn at Q = is a finite-size effect 
which also occurs for the pure spin model and rapidly de- 
creases with increasing N [^o|; for N — > oo, EJ(tt) and 
Eq (0) become degenerate. 

The Peierls ordering structure of the groundstate can 
be demonstrated in two complementary ways. The mag- 
netic ordering pattern of an alternating local singlet for- 
mation is captured by the order parameter 



^ = (^E(- 1 ) i ( s -i- s '- s ^ s ^)) 



(4) 



D is the proper choice for the Heisenberg spin chain 
model (|l|). For the translationally invariant quantum 
spin-phonon Hamiltonian ^ the Peierls distortion of the 
lattice is reflected in the spatial structure of the displace- 
ment correlation function 



C hl = ((b 1 +b+) (fei + ftt))-^ 



(5) 




FIG. 4. Displacement correlation function Ci,i as defined 
in the text for different phonon frequencies and spin lattice 
coupling constants for a chain with TV = 12 sites. 



and long bonds and thus alternating strong and weak 
AF exchange interactions, i.e. a dimerized groundstate. 
Naturally the dimerization is enhanced (weakened) by in- 
creasing the spin-phonon coupling (phonon frequency). 
We note that for any finite phonon frequency there exists 
a critical coupling constant g c plj ; in the infinite sys- 
tem quantum lattice fluctuations will destroy the Peierls 
dimerization for couplings g < g c . For non-adiabatic 
phonon frequencies hu; ~ J also the critical coupling g c 
is comparable to J - a situation which matches the strong 
magnetoelastic coupling and the SP physics in CGO. 

In order to quantify these results we obtain the magni- 
tude of the Peierls dimerization order parameter S from 
the static (lattice) structure factor at wavenumber q = n 



52 = I^E^i^-^ 

id 



(6) 



The alternatingstructure of the correlation function Ci,, 
shown in Fig. implies the Peierls formation of short 



where Ui = bi + . Like in the ordinary Peierls phe- 
nomenon, a finite dimerization 5 > necessarily leads 
to a gap A ST in the magnetic excitation spectrum. For 
evaluating the relation between the dimerization and the 
resulting magnitude of the spin triplet excitation gap 
we keep the phonon frequency fixed, vary the coupling 
strength g and calculate for each parameter set A ST (see 
Fig. |E[), and the dimerization S from Eq. (^). 

The results for the static (spin-only) and the dynamic 
spin-phonon model are compared in Fig. [5] for three dif- 
ferent chain lengths. For vanishing lattice dimerization, 
i.e. in the absence of any spin-phonon coupling [g = 0), 
the results for the static and the dynamic model naturally 
agree (note that for (5 = there remains a spin excitation 
gap even for N — > 00 due to the frustration driven sin- 
glet dimer ordering D ^ 0). The important message of 
Fig. H is that a substantial enhancement of the dimeriza- 
tion is needed to achieve a given triplet excitation gap, if 
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performed on the vector-parallel supercomputer Fujitsu 
VPP700 at the LRZ Munchen. 
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FIG. 5. Triplet excitation gap A ST versus dimerization 
5 for the static spin-only (uj — 0, o □ o) and the dynamic 
spin-phonon model (full symbols) with phonon frequency 
uj — 1. for different chain lengths at a = 0.36. The inset shows 
the results for the dynamic eight-site model (filled symbols; 
uj = 0.1 [0.3] triangles up [down], 1.0 (•)) in comparison to 
the static case (o). 



the dynamic phonons have frequencies comparable to J. 
This partially resolves the quantitative problem encoun- 
tered within the alternating and frustrated Heisenberg 
spin chain model (|l|) , because fixing parameters J and a 
by fitting the high temperature susceptibility and 6 by 
matching A ST to the INS data in CGO, the dimeriza- 
tion par ameter 6 is underestimated by roughly a factor 
of 5 [|19|,j32|| . The inset in Fig. || shows the evolution of 
the A 6J vs. S results with increasing phonon frequency 
to at a fixed chain length. With decreasing uj the data 
smoothly merge with the result for the static spin-only 
model. 

In conclusion, we have numerically analyzed the Peierls 
dimerized groundstate and the spin triplet excitations 
in a Heisenberg chain model with dynamic spin-phonon 
coupling. The magnetic excitations inherently include a 
local lattice distortion requiring a multi-phonon mode 
treatment of the lattice degrees of freedom. A spin- 
phonon coupling induced flattening of the triplet dis- 
persion is avoided for a strong spin-phonon coupling in 
the non-adiabatic regime. The non-adiabatic phonon 
dynamics strongly renormalizes the magnetic excitation 
spectrum and the dimerization dependence of the triplet 
excitation gap. These features are in accordance with 
experimental data on CGO suggesting the necessity for a 
non-adiabatic spin-phonon approach to the SP physics 
in this material. 
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the SP phase of CGO. The numerical calculations were 
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